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Our universe contains a great number of extremely compact and massive objects which are gen- 
erally accepted to be black holes. Precise observations of orbital motion near candidate black holes 
have the potential to determine if they have the spacetime structure that general relativity demands. 
As a means of formulating measurements to test the black hole nature of these objects, Collins and 
Hughes introduced "bumpy black holes": objects that are almost, but not quite, general relativity's 
black holes. The spacetimes of these objects have multipoles that deviate slightly from the black 
hole solution, reducing to black holes when the deviation is zero. In this paper, we extend this work 
| in two ways. First, we show how to introduce bumps which are smoother and lead to better behaved 

, orbits than those in the original presentation. Second, we show how to make bumpy Kerr black 

holes — objects which reduce to the Kerr solution when the deviation goes to zero. This greatly 
' extends the astrophysical applicability of bumpy black holes. Using Hamilton- Jacobi techniques, we 

show how a spacetime's bumps are imprinted on orbital frequencies, and thus can be determined by 
■ measurements which coherently track the orbital phase of a small orbiting body. We find that in the 

Qweak field, orbits of bumpy black holes are modified exactly as expected from a Newtonian analysis 
of a body with a prescribed multipolar structure, reproducing well-known results from the celestial 
£ — ■ mechanics literature. The impact of bumps on strong-field orbits is many times greater than would 

be predicted from a Newtonian analysis, suggesting that this framework will allow observations 
to set robust limits on the extent to which a spacetime's multipoles deviate from the black hole 
expectation. 



PACS numbers: 04.25.Nx, 04.30.Db, 04.70.Bw 



I. INTRODUCTION 



(N 

, A. Motivation: Precision tests of the black hole hypothesis 

£ — , . Though observations of gravitational systems agree well with the predictions of general relativity (GR), the most 
detailed and quantitative tests have so far been done in the weak field. ("Weak field" means that the dimensionless 
Newtonian potential = GM/rc 2 <C 1, where M is a characteristic mass scale and r a characteristic distance.) 
■ This is largely because many of the most precise tests are done in our solar system (e.g., [l|). Even the celebrated 
Q\ [ tests which use binary neutron stars (e.g., (H) are essentially weak-field: for those systems, M ~ severalM©, r = 
. orbital separation ~ several i?©, so <f> ~ a f ew x GMq/Rqc 2 ~ a few x 10 -6 . 

This situation is on the verge of changing. Observational technology is taking us to a regime where we either are 
or soon will be probing motion in strong gravity, with <fi > 0.1. Examples of measurements being made now include 
radio studies of accretion flows near the putative black hole in our galactic center [U and x-ray studies that allow 
precise measurements of accretion disk geometries [p, Hj]. Future measurements include the possible discovery of a 
black hole-pulsar system, perhaps with the Square Kilometer Array Q , and gravitational- wave (GW) observations of 
small bodies spiraling into massive black holes due to the backreaction of GW emission 0]. 

For weak-field studies, a well-developed paradigm for testing gravity has been developed. The parameterized post- 
Newtonian (PPN) expansion [H, Q quantifies various measurable aspects of relativistic gravity. For example, the PPN 
parameter 7 (whose value is 1 in GR) quantifies the amount of spatial curvature produced by a unit of rest mass. 
Other PPN parameters quantify a theory's nonlinearity, the degree to which it incorporates preferred frames, and 
the possible violation of conservation laws. See Ref. 0, Chap. 4 for a detailed discussion. Unfortunately, no similar 
framework exists for strong-field studies. If we hope to use observations as tools for testing the nature of strong 
gravity objects and strong-field gravity, we need to rectify this. 

Black holes are of particular interest for studying strong-field gravity. Aside from having the strongest accessible 
gravitational fields of any object in the universe 1 , within GR they have an amazingly simple spacetime structure: the 



More accurately, they have the largest potential <f> ~ GM/ Rc 2 . One may also categorize weak or strong gravity using spacetime curvature 
or tides; arguably this is a more fundamental measure for assessing whether GR is likely to be accurate or not. From this perspective, 
black holes are actually not such "strong gravity" objects; indeed, the tidal field just outside a 10 s Mq black hole's event horizon is not 
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"no-hair" theorems [ljj, HH, [H, [H, (3 guarantee that the exterior spacetime of any black hole is completely described 
by only two numbers, its mass M and spin parameter a. Any deviation from that simplicity points to a failure either 
in our understanding of gravity or in the nature of ultracompact objects. 

Recent work by Brink [l5| has reviewed in detail the challenges involved in testing the spacetime structure of 
massive compact bodies when one relaxes the assumption that the spacetime is Kerr. As in Brink's discussion, we 
will focus on spacetimes that are stationary, axisymmetric, and vacuum in some exterior region. An especially useful 
parameterization is provided by the work of Geroch [l6j and Hansen [l7j . They demonstrate that such a spacetime 
is completely specified by a set of mass multipole moments M; and current multipole moments Si. For a fluid body, 
Mi describes the angular distribution of the body's mass, and Si describes the angular distribution of mass flow. In 
Newtonian theory, M; labels a piece of the gravitational potential that varies as Yi (9, cf>)/r l+1 . There is no Newtonian 
analog to Si . In the weak- field limit of GR, they label a magnetic-like contribution to gravity and are reminiscent of 
magnetic moments. 

In general, the moments {Mi, Si} can be arbitrary. For black holes, they take special values: a Kerr spacetime has 

m 

Mi+iSi = M(ia) 1 , (1.1) 

(in units with G = 1 = c; we neglect the astrophysically uninteresting possibility of the hole having macroscopic 
charge). In other words, only two moments out of the set {Mi, Si} are independent: Mo = M, and Si — aM. Once 
those two moments have been specified, all other moments are fixed by Eq. (jl.ip if the spacetime is Kerr. 

Testing the hypothesis that an object is a Kerr black hole can thus be framed as a null experiment. First, measure 
the putative black hole spacetime's multipoles. Using the moments Mq and Si, determine the parameters M and a. 
If the spacetime is Kerr, all moments for I > 2 must be given by Eq. The null hypothesis is that any deviation 

from those Kerr moments is zero. Failure of the null hypothesis means the black hole candidate is not a Kerr black 
hole, and may indicate a failure of strong-field GR. 

B. Bumpy black holes: Past work 

Performing this null experiment requires strong-field spacetimes with multipoles differing from those of black holes. 
Given the tremendous success of general relativity and the black hole hypothesis at explaining a wide span of data, 
a reasonable starting point is to consider spacetimes whose multipoles deviate only slightly. Imagine starting with a 
spacetime whose Geroch-Hansen moments satisfy 

Mi + iSi = M{ia) 1 + 5Mi + iSSi . (1.2) 

The null hypothesis — black hole candidates are described by GR's Kerr metric — means 5 Mi = 0, 5 Si — for all I. 

To this end, Collins and Hughes [l8[ (hereafter CH04) introduced the bumpy black hole: a spacetime that deviates 
in a small, controllable manner from the exact black holes of GR. By construction, the bumpy black hole includes 
"normal" black holes as a limit. This is central to testing the black hole hypothesis by a null experiment. Spacetimes 
of other proposed massive, compact objects (for example, boson stars 0, H(|) typically do not include black holes 
as a limiting case. This limits their utility if black hole candidates are, in fact, GR's black holes. Measurements of 
black holes using observables formulated in a bumpy black hole spacetime should simply measure the spacetime's 
"bumpiness" (defined more precisely below) to be zero. 

Though a useful starting point, the bumpy black holes developed in CH04 had two major shortcomings. First, 
the changes to the spacetime that were introduced to modify its multipole moments were not smooth. The worked 
example presented in CH04 is interpreted as a Schwarzschild black hole perturbed by an infinitesimally thin ring of 
positive mass around its equator and by a pair of negative mass infinitesimal points near its poles. Though this changes 
the spacetime's quadrupole moment (the desired outcome of this construction), it gives the spacetime a pathological 
strong-field structure. This is reflected in the fact that non-equatorial strong-field orbits are ill-behaved in the CH04 
construction [2lj . 

Second, CH04 only examined bumpy Schwarzschild black holes. Although their approach differed in many details 
from that used in CH04, Glampedakis and Babak [22| (hereafter GB06) rectified this deficiency with their introduction 
of a "quasi-Kerr" spacetime. Their construction uses the exterior Hartle-Thorne metric [23], [2J| describing the exterior 
of any slowly rotating, axisymmetric, stationary body. It includes the Kerr metric to 0(a 2 ) as a special case. 



much different from the tidal field at the surface of the Earth. See Ref. |f| for further discussion of this point. 



3 



Identifying the influence of the quadrupole moment in the Hartle-Thorne form of the Kerr metric, they then introduce 
a modification to the full Kerr spacetime that changes the black hole's quadrupole moment from its canonical value. 

C. This analysis 

The goal of this work is to rectify the deficiencies of CH04 and to extend the bumpy black hole concept to spinning 
black holes. Dealing with the non-smooth nature of the bumps presented in CH04 is, as we show in Sees. HJ IIV1 and 
IVl quite straightforward. It simply requires introducing perturbations to the black hole background that are smooth 
rather than discontinuous. In essence, rather than having bumps that correspond to infinitesimal points and rings, 
the bumps we use here are smeared into pure multipoles. 

Extending CH04 to spinning black holes is more of a challenge. Given that GB06 already introduced a bumpy- 
black-hole-like spacetime that encompasses spacetimes with angular momentum, one might wonder why another 
construction is needed. A key motivation is that we would like to be able to make an arbitrary modification to a black 
hole's moments. Although showing that a black hole candidate has a non-Kerr value for the quadrupole moment 
would be sufficient to falsify its black hole nature (at least within the framework of GR), one can imagine scenarios 
in which the first L moments of a black hole candidate agree with the Kerr value, but things differ for I > L. For 
example, Yunes et al. have shown that in Chcrn-Simons modifications to GR, slowly rotating black hole solutions 
have the multipolar structure of Kerr for I < 4, but differ for I > 4 [25^ . There are many ways in which black hole 
candidates might differ from the black holes of general relativity; we need to develop a toolkit sufficiently robust that 
it can encompass these many potential points of departure. 

Our technique for making bumpy Kerr holes is based on the Newman- Janis algorithm [2q |. This algorithm trans- 
forms a Schwarzschild spacetime into Kerr by "rotating" the spacetime in a complex configuration space 2 . In this 
paper, we construct bumpy Kerr black holes by applying the Newman- Janis algorithm to bumpy Schwarzschild black 
holes. The outcome of this procedure is a spacetime whose mass moments are deformed relative to Kerr. When 
the bumpiness is set equal to zero, we recover the Kerr metric. A companion paper (28j examines the multipolar 
structure of bumpy Kerr black holes in more detail, demonstrating that one can construct a spacetime in which the 
mass moments deviate arbitrarily from Kerr (provided that they are small) . 

Once one has constructed a bumpy black hole spacetime, one then needs to show how its bumps are encoded in 
observables. The most detailed quantitative tests will come from orbits near black hole candidates. As such, it is 
critical to know how orbital frequencies change as a function of a spacetime's bumpiness. More generally, we are 
faced with the problem of understanding motion in general stationary axisymmetric vacuum spacetimes. Brink (29| 
has recently published a very detailed analysis of this problem, with a focus on understanding whether and for which 
situations the spacetimes admit integrable motion. She has found evidence that geodesic motion in such spacetimes 
may, in many cases, be integrable. If so, the problem of mapping general spacetimes (not just "nearly black hole" 
spacetimes) may be tractable. Gair, Li, and Mandel [3fj| have similarly examined orbital characteristics in the Manko- 
Novikov spacetime [3~lj . which has a particular tunable non-Kerr structure. They show how orbits change in such 
spacetimes, and how its bumpiness colors observable characteristics. 

For this analysis, we confine ourselves to the simpler problem of motion in bumpy black hole spacetimes, addressing 
this challenge in Sees. IIVI and [V] using canonical perturbation theory. As is now well known (and was shown rather 
spectacularly by Schmidt [32j), Hamilton- Jacobi methods let us write down closed-form expressions for the three 
orbital frequencies (f2 r , fi , fi*) which completely characterize the behavior of bound Kerr black hole orbits. Since 
bumpy black hole spacetimes differ only perturbatively from black hole spacetimes, canonical perturbation theory 
lets us characterize how a spacetime's bumps shift those frequencies, and thus are encoded in observables. Similar 
techniques were used in GB06 to see how frequencies are shifted in a quasi-Kerr metric, and were also used by Hinderer 
and Flanagan [33] in a two-timescale analysis of inspiral into Kerr black holes. 

D. Organization and overview 

We begin in Sec. [Til with an overview of the spacetimes that we study. We start with the axially symmetric 
and stationary Weyl line element. We first review the Einstein field equations in this representation, introduce the 
Schwarzschild limit, and then describe first order perturbations. The spacetime's bumpiness is set by choosing a 



2 This complex "rotation" is made most clear by framing the discussion using the Ernst potential |27| . Transforming from Schwarzschild 
to Kerr corresponds to adding an imaginary part to a particular potential. Further discussion is given in a companion paper [28l] , 
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function ipi which controls how the spacetime deviates from the black hole limit. We initially leave ipi arbitrary, 
except for the requirement that it be small enough that terms of order (V>i) 2 can be neglected. Later in the paper, 
we will take ?pi to be a pure multipole in the Weyl representation. Following our discussion of bumpy Schwarzschild 
spacetimes, we show how to use the Newman- Janis algorithm to build bumpy Kerr black holes. 

We study geodesic motion in these spacetimes in Sec. IIII1 We begin by reviewing the most important properties of 
normal black hole orbits, reviewing Kerr geodesies in Scc. lIII A"l and then describing how to compute orbital frequencies 
using Hamilton- Jacobi methods in Sec. lIII 51 The discussion of frequencies is largely a synopsis of Schmidt's pioneering 
study, Ref. [321 ] . We then show how these techniques can be adapted, using canonical perturbation theory, to bumpy 
black holes fSec. Hn~C|) . Canonical perturbation theory requires averaging a bump's shift to an orbit's Hamiltonian. 
This averaging was developed in Ref. (4p| , and is summarized in Appendix [A1 

We present our results for bumpy Schwarzschild and bumpy Kerr orbits in Sees. IIVI and fVl respectively. We take 
tpi to be a pure multipole in the Weyl sector, construct the spacetime, and numerically compute the shifts in SI 1 ", f] e , 
and SI*. Detailed results are given for 1 = 2,1 = 3, and / = 4. There is no reason in principle to stop there, though 
the results quickly become repetitive. Section |IV] gives our results for bumpy Schwarzschild black holes, and Sec. \V\ 
gives results for bumpy Kerr. 

A few results are worth highlighting. First, we find that the exact numerical results for frequency shifts correctly 
reproduce the Newtonian limit as our orbits are taken into the weak field. (We develop this limit in some detail in 
Appendix [Bj to facilitate the comparison.) The frequency shifts are substantially enhanced in the strong field. The 
shift to the radial frequency is particularly interesting: it tends to oscillate, shifting between an enhancement and a 
decrement as orbits move into the strong field. This behavior appears to be a robust signature of non-Kerr multipole 
structure in black hole strong fields. 

Interestingly, it turns out that black hole spin does not have a very strong impact on the bumpiness-induced shifts 
to orbital frequencies. Spin's main effect is to change the location of the last stable orbit. For large spin, orbits reach 
deeper into the strong field, amplifying the bumps' impact on orbital frequencies. Aside from the change to the last 
stable orbits, the impact of a particular multipolar bump looks largely the same across all spin values. Some examples 
of the frequency shifts we find are shown in Figs. [T] (/ = 2, Schwarzschild), =4, Schwarzschild) and [3] (I = 2, 
Kerr). (We have no figures for / = 3 since the secular shifts to orbital frequencies are zero in that case, as they are for 
all odd values of I. We also do not show results for I = 4 bumps of Kerr black holes since, as discussed above, they 
are very similar to the I = 4 Schwarzschild results.) 

Finally, we summarize our analysis and suggest some directions for future work in Sec. I VII Among the points we 
note are that, in this analysis, we only consider changes to the mass moments of the black hole spacetime. Adding 
"current-type" bumps to a spacetime is discussed in a companion paper *28j|. We also do not discuss in this analysis 
the issue of measurability. Turning these foundations for mapping the multipole moment structure of black holes into 
a practical measurement program (for instance, via gravitational-wave measurements, timing of a black hole-pulsar 
binary, or precision mapping in radio or x-rays of accretion flows) will take a substantial effort. 

Throughout this paper, we work in geometrized units with G = c = 1; a useful conversion factor in these units 
is 1 M© = 4.92 x 10 -6 seconds. When we discuss bumpy black holes, we will always use a "hat" accent to denote 
quantities which are calculated in the pure black hole background spacetimes. For example, an orbital frequency is 
written SI = fl + <5S7; SI is the frequency of an orbit in the black hole background, and <5S! denotes the shift due to the 
black hole's bumps. 



II. BLACK HOLE AND BUMPY BLACK HOLE SPACETIMES 



We begin with general considerations on the spacetimes we consider. Since we focus on stationary, axisymmetric 
spacetimes, the Weyl metric [36| is a good starting point: 



ds 2 = -e 2 Ut 2 + e 2 ^(dp 2 + dz 2 ) + e- 2 ^p 2 d<f 
The nontrivial vacuum Einstein equations for this metric are given by 





dj 



d 2 ip 1 dtp d 2 ip 
dp 2 p dp dz 2 



dtp 
dp 



dtp 
dz 



df ^ 

dz dp dz 



(2.1) 

(2.2) 
(2.3) 
(2.4) 
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Equations (|2.2[) - (|2.4p will be our main tools for building bumpy black hole spacetimes. We will put -0 = ipo + ip\ , 
7 = 7o + 7ii w itb tpi/ipo <C 1, and 71/70 -C 1. Before specializing to black hole backgrounds, note that Eq. (|2.2[) is 
simply Laplace's equation. The functions -01 can thus very conveniently be taken to be harmonic functions. This is 
key to smoothing out the spacctime's bumps and curing one of the deficiencies of CH04. 



A. Schwarzschild and bumpy Schwarzschild 

We begin by building a bumpy Schwarzschild black hole. The Schwarzschild metric is recovered from Eq. (|2.ip 
when ipi = 71 = and we set 

0o = lntanh(u/2) , (2.5) 
1 , / sin 2 v 



7o = --m 1 + —^ . (2.6) 
* \ sinn it J 

The prolate spheroidal coordinates (u,v) are a remapping of the coordinates (p, z) used in Eq. (|2.1|) : 

p = M sinh u sin v , (2.7) 
z = M cosh it cos v . (2-8) 

(Below we will remap these to the familiar Schwarzschild form.) Expanding the Einstein equations Eqs. (|2.2p - (12. 4|) 
about these Schwarzschild values to leading order in 0i and 71, the perturbations must satisfy 

VVi = 0, (2.9) 
371 2[ta,nv(dipi / du) + ta,nhu(dipi/dv)] 



du sinh u(coth it tan v + tanhu cot v) 

dji 2[t&nv(dijji / dv) — tanhu(dipi/du)] 



(2.10) 
(2.11) 



dv sinh u(coth it tan v + tanhu cot v) 

As discussed in CH04, Eqs. (|2.10[) and (|2.11[) actually overdetermine the solution; we will use Eq. (|2.11[) to calculate 
71. [Note also that the tanv in the numerator of Eq. (|2. 10|1 is incorrectly written cotv in CH04.] 

To connect this to the variables we will use later in the paper (and to put it in a more familiar form), we make a 
final change of coordinates, putting 

r = 2Af cosh 2 (it/2) , (2.12) 

6 = v , (2.13) 

so that 



p = rsin6>Wl - — , (2.14) 
z = (r-M)cos9 . (2.15) 



The spacetime then becomes 



ds 2 = -e 2 ^ 1^1 ~ ™\ dt 2 + e 2 ^- 2 ^ 1 - dr 2 + r 2 e 2 ^~ 2 ^de 2 + r 2 sin 2 9 e - 2 ^d(b 2 

= (g a p + b aP )dx a d X P . (2.16) 

Although we have left the potentials ip\ and 71 in exponential form, these quantities must be expanded to first 
order, since we solve for them using linearized Einstein equations. We use the exponential form only for notational 
convenience. On the second line, g a p is the Schwarzschild metric, and 

ht = -20! - , (2.17) 

b rr = (271-2^1) f 1 - — ) . ( 2 - 18 ) 

b ee - (2 7l -20!)r 2 , (2.19) 
6 00 = ~20ir 2 sin 2 6» . (2.20) 

All other components of b a p are zero. We clearly recover the normal Schwarzschild black hole when ip\ — > 0, 71 — > 0. 
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B. Kerr and bumpy Kerr 

We now use the Newman- Janis algorithm [2t| to transform bumpy Schwarzschild into bumpy Kerr. We begin with 
the bumpy Schwarzschild metric written in prolate spheroidal coordinates: 

ds 2 = -e 2 ^ 1 tanh 2 (w/2)fft 2 + e 27l - 2 *4M 2 cosh 4 (u/2) (du 2 + dv 2 ) + e" 2 *4M 2 cosh 4 (w/2) sin 2 v d(/) 2 . (2.21) 

The first step in the Newman-Janis algorithm uses the fact that the metric can be written in terms of a complex null 
tetrad with legs l^, n v , m u : 

gT = -Fn u - Fn^ + m^m u + m u rh^ ; (2.22) 
an overbar denotes complex conjugate. The legs are given by 

l" = e^' 1 coth 2 (u/2) 6£ + -jjU^-^cschu^ , (2.23) 

= -e-^S? e^-^cschu tanh 2 (u/2) #f , (2.24) 

2 * 2M v ' ' " 

m" = 2 ^L^ e^sech 2 (M/2) (e~ 71 5ff + zcscu<5g) . (2.25) 

We use the Kronecker delta 5^ to indicate components. Next follows the key step of the Newman-Janis algorithm: 
We allow the coordinate u to be complex, and rewrite l^, n^, and m M as 

= e - * 1 5 * 5 +e* i -~* — , (2.26) 

tanh (u/2) + tanh (u/2) Mvcoshu coshw — 1 

1 , „ , [tanh 2 (u/2) + tanh 2 (u/2)] 

n » = -4:1 g? _ e *i-7i i \ 7 7 ^=4- 6£ , (2.27) 

2 4Mv cosh u cosh it — 1 

= —i—e^sech 2 ^) (e -71 ^ + icscu^) . (2.28) 

Notice that we recover the original tetrad when we force u = u. Further discussion of this seemingly ad hoc procedure 
(and an explanation of how it uniquely generates the Kerr spacetime) is given in Ref. [37] . 
Next, change coordinates: Rewrite the tetrad using (U, r, 9, cf>), given by 

U = t- 2M cosh 2 (u/2) - 2M In [sinh 2 (it/2)] - iacos9 , (2.29) 
r = 2M cosh 2 (u/2) + ia cos 9, (2.30) 
9 = v . (2.31) 

The axial coordinate <f> is the same in both coordinate systems. At this point, a is just a parameter. The result of 
this transformation is 



( 2Mr\ _1 

e^-^) f 1- — — J <y^ + e* 1- ^ , (2.32) 



I( e -*+ e *-»)#-±e*-» (l - ^) 



5? , (2.33) 



1a sin ( 



V / 2('" + cos 

where £ = r 2 + a 2 cos 2 0. Making one further coordinate transformation 



((5£-^) + <^ + e 71 icsc6l<5£ , (2.34) 



dU = dt- r ~T Q dr , dcf> = #' - ^dr , (2.35) 
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gives us a bumpy Kerr black hole metric in Boyer-Lindquist coordinates: 

2Mr 



ds 1 



=,2^1 



1 



^ + e^(l-e^) 4a2Mr „ Sin2g dtdr- 



2Mr\ 



1 + e~ 271 (l - 2e 71 ) 



AS 
a 2 sin 2 9 



,4^1-471 _ e 7iy 



,4a 4 M 2 r 2 sin 4 ^" 1 
A 2 E 2 



dr 2 



-2(1 - e 7l )asin 2 1 
_ e 2 7l -2^ 1E(i0 2 + A 



-2i/>i 



1 - 



2Mr 



-2i/>i 



s / 

2Mr\ 



\ _ e 2Vi-2 7 i 



^i-27! 



4a 2 M 2 r 2 sin 2 
AE(E - 2Mr) 

4a 2 M 2 r 2 sin 2 



AE(E - 2Mr) 



(2.36) 



where A = r 2 — 2Mr + a 2 , and we have dropped the prime on <fr. Notice that Eq. (|2.36p reduces to the "normal" Kerr 
black hole metric when ipi — > 0, 71 — > 0; the parameter a is seen to be the specific spin angular momentum, \S\/M. 
A companion paper (28j examines the multipoles of this spacetime for particular choices of ipi and demonstrates that 
it corresponds to Kerr with some moments set to the "wrong" values. The Newman- Janis algorithm applied to the 
bumpy Schwarzschild black hole produces a bumpy Kerr black hole. 

Writing the bumpy Kerr metric in the form g a @ = g a p + b a p, we read out of Eq. (|2.36p 



hr 

bt</> 

b r <p 
bee 

bthtb 



-2 1 



2Mr\ 



S 

2a 2 Mr sin 2 9 



4>i 



-71- 
(71 - 
2(71 

71 
2(71 
(71 



2w 



AS 

% 2aMrsin 2 6» 



^1) 



1 



2Mr 



Aa 2 M 2 r 2 sin 2 9 
AS(S - 2Mr) 



E 

Vi)s 

. 8a 2 M 2 r 2 sin 2 9 
W AE(E - 2A/r) ~ 2 ^ 



a sin 2 



2Mr\ 1 



A sin 2 6> 



(2.37) 
(2.38) 
(2.39) 
(2.40) 

(2.41) 
(2.42) 
(2.43) 



Other components are related by symmetry or zero. By inspection, we can see that b a p — ► as ifii — * 0, 71 — * 0. 

Before moving on, we summarize. To build a bumpy black hole spacetime, we first select a function %p\ which 
satisfies the Laplace equation (|2.2[) . We find the function 71 which satisfies Eqs. (I2.3P and (|2.4p . and then apply the 
Newman- Janis algorithm to "rotate" the spacetime to non-zero a. The result is given by Eqs. (|2.36p - (|2.43D . 



III. MOTION IN BLACK HOLE AND BUMPY BLACK HOLE SPACETIMES 



We now discuss motion in these spacetimes. Our focus will be computing the frequencies associated with oscillations 
in the radial coordinate r, the polar angle 9, and rotations in <f> about the symmetry axis. These frequencies are 
typically the direct observables of black hole orbits; it is from measuring these frequencies (or the evolution of these 
frequencies if the orbit evolves) that one can hope to constrain the properties of black hole candidates. 



A. Geodesies of Kerr black holes 



As background to our discussion of black hole orbital frequencies, we first briefly review the equations governing 
black hole geodesic orbits, and some useful reparameterizations for practical computations. One goal of this discussion 
is to introduce certain quantities which we will use in the remainder of this paper. 

As first recognized by Carter [38[, geodesic motion for a test mass m in a black hole spacetime is separable with 
respect to Boyer-Lindquist coordinates t, r, 9, and <j). The test body's motion is then completely described by four 
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first-order ordinary differential equations: 

2 



2v^2 
m 2j 



2v^2 

m 2j 



[(r 2 +a 2 )£-aL 2 ] - A [m 2 r 2 + (L z 



aEf + Q] 



mS — = 



R(r) , 

Q - cot 2 9L 2 Z - a 2 cos 2 
0(9), 

esc 2 6L z +aE 
$0,0) , 



£ 2 ) 



r 2 + a 2 



-1 - 



(r 2 + a 2 ) 2 



a 2 sin 2 i 



a 2 £ 2 
A 



■ aL, 1 - 



T(r,0) 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



As in the previous section, A = r 2 — 2Mr + a 2 and £ = r 2 + a 2 cos 2 9; r is proper time along the test body's 
worldline. In developing these equations, one isolates four constants of the motion. One is simply the rest mass itself, 
m 2 = —p^p p ; this motivates the definition of the Hamiltonian for test body motion, 



(3.5) 



where the 4-momentum components p^ = mg lLV [dx v j dr). The other constants are the energy E, axial angular 
momentum L z , and "Carter constant" Q, given by 



E = -p t 



Q 



cos 2 9 [a 2 (m 2 -E 2 )+ esc 2 9L 2 Z ] 



(3.6) 
(3.7) 
(3.8) 



Given a choice of the constants (E,L Z ,Q) and a set of initial conditions, Eqs. (|3.1[) - p.4| completely describe the 
geodesic motion of a test body near a Kerr black hole. The equations for r and 9 can present some problems, however, 
since their motion includes turning points where dr/dr and d9/dr pass through zero and switch sign. To account for 
this behavior, it is convenient to reparameterize these motions using angles ip r (for the radial motion) and \ (for the 
polar motion) which smoothly vary from to 2ir as the motion oscillates between its extremes. 
Consider first the radial motion. We define 



pM 



1 + e cos i/v 



(3.9) 



The constants p and e are the orbit's semi-latus rectum and eccentricity, respectively. Substituting into Eq. (|3.ip . it is 
simple to develop an equation governing ip r . Periapsis and apoapsis are given by r p — pMj (1 + e) and r a = pM / (1 — e), 
respectively. This allows us to relate the constants p and e to the constants E, L z , and Q: They are the outermost 
radii for which the radial "potential" goes to zero, 



For the polar motion, we note that Eq. 



R(r p ) = R{r a ) = 
can be written 



, 'd9\ 2 _ z 2 [a 2 {m 2 -E 2 )} - z [Q + L 2 Z + a 2 (m 2 - E 2 )} + Q 



(3.10) 



(3.11) 



where we have introduced z = cos 2 9. Denote by z± the two roots of the quadratic on the right-hand side of Eq. 
(|3.1ip . Turning points of the 9 motion correspond to z = z_ , the smaller of these roots. (The root z + is greater than 
1, and does not correspond to a turning point.) Transforming back to the angle 9, we find that the minimum polar 
angle reached by the orbit is 



cos 6> min = ^Jz~Z . 



(3.12) 
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(The maximum angle is # m ax = it — O m - ln .) A useful reparameterization of the 8 coordinate is 

COS = COS 6> m in cos % \ (3.13) 

where x accumulates like the angle tp r . By substitution in Eq. p.2p , we can easily develop an equation governing the 
evolution of \. 

Before moving on, we note that Eqs. (|3.10|) and (|3.12|) allow us to map from the parameters (E, L z , Q) to (p, e, 6 m i n ). 
[Schmidt [32J in fact gives an analytic solution for (E,L Z ,Q) as functions of [p, e, Q m in)-\ Up to initial conditions, 
either parameterization thus completely specifies an orbit. We will flip between these parameterizations as convenient. 



B. Orbital frequencies for black holes 

To frame our discussion, we begin by examining the frequencies of motion for "normal" black hole spacetimes. Our 
discussion here closely follows that given by Schmidt [12] , which uses Hamilton- Jacobi methods to compute black hole 
orbital frequencies. A useful starting point is to note that in separating the coordinate-space motion, one not only 
identifies constants of the motion, but also the action variables 



J r ee -L Ip r dr=- r ^1 dr , (3.14) 



2n J 7T 

i r ,„ 2 



tt/2 



jg ee —j Pg d0 = - I ^/eiejdo, (3.15) 



1 

27 



(3.16) 



It is also useful to define 



J t = -E. (3.17) 

This is a slight abuse of the notation since geodesic motion is not cyclic in t (and hence we cannot define Jt as a 
closed integral over time), but is convenient for reasons we will illustrate shortly. 

At least formally, we can now reparameterize our Hamiltonian (|3.5[) in terms of the action variables J M . Let us write 
the Hamiltonian so reparameterized as 7i^ aa \ By Hamilton- Jacobi theory, the orbital frequencies are the derivatives 
of 7i( aa ) with respect to the action variables: 

<9H (aa) 

W = _. (3.18) 

For black hole orbits, we cannot explicitly reparameterize the Hamiltonian in this way, making it difficult to calculate 
the orbital frequencies. However, by using the chain rule, we can re-express (13.181) in terms of derivatives that are not 
so difficult to explicitly write out. Following Ref. 32] (modifying its notation slightly), we put 

P = (H,E,L Z ,Q) . (3.19) 

Define the matrices A and B, whose components are 

p _ dP„ p _ dJa , . 

A a - dJ ^ , B a - . (3.2U) 

By the chain rule, these matrices have an inverse relationship: 

AjBfP = S a t . (3.21) 

The components of the matrix A are directly related to the frequencies we wish to compute. In particular, since Pt 
is just the invariant Hamiltonian, muj 1 = dPt/dJi = At 1 . However, the components of the matrix B are written in a 

way that is fairly easy to work out. We exploit this to write mu' = (B -1 )^ , from which we find 

r dJe/dQ 



{dJ r /dH){dJ e /dQ) - {dJ r /dQ)(dJ e /dH) ' 

-dJ r /dQ 

{dJ r /dH){dJ /dQ) - {dJ r /dQ)(dJ e /dH) ' 



(3.22) 
(3.23) 



4 = (dJ r /dQ)(dJ e /dL z ) - (dJ r /dL z )(dJ e /dQ) 
mU} {dJ r /dH){dJg/dQ)-{dJ r /dQ){dJ e /dH) ' 1 ' ' 
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The partial derivatives of J r and Jg appearing here are given by 

„2 



d,J r 
BH 
dJ r 
~dQ 

dJ r 
dL~ 

dJ r 
~dE 

dU 

dQ 

dJ e 
dL z 

dJe_ 

dE 



1 
2^ 

1 

47T 
1 

1 

2^ 
1 

2^ 
1 

4tt 

1 

1 

27 



1 



r{rL z - 2M(L Z - aE)) 



dr , 



r [rE(r 2 + a 2 ) - 2Ma(L z - aE)] 



dr 



a 2 cos 2 I 



1 



■d6> 



-,de 



VW) 

L z cot 2 6> 



VW) 
a 2 E cos 2 6* 



dff 

de . 



(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 



(The derivatives dJ r .e/dE will be needed for a quantity we introduce below.) Schmidt [32j combines these results to 
give closed-form expressions for the three frequencies uf' we will not repeat these expressions here. 

The frequencies co r ' e '^ are conjugate to the orbit's proper time; they would be measured by an observer who rides 
on the orbit itself. For our purposes, it will be more useful to convert to frequencies conjugate to the Boyer-Lindquist 
coordinate time, describing measurements made by a distant observer. The quantity 3 



r 



performs this conversion; the frequencies fl r ' 
and (|3~2Tj) . we find mT = (B~ 



1 cW (aa) 1 <9H (aa) 



(3.33) 



t , or 



to dJt to dE 

,r,e,<p jY are f observational relevance. Going back to Eqs. (|3.20|) 



toT = 



(dJ r /dE){dJg/dQ) - (dJ r /dQ){dJg/dE) 

(dJr/dH)(dJ e /dQ) - (dJ r /dQ)(dJe/dH) 



(3.34) 



In the discussion that follows, it will be useful to have weak- field (j> 3> M) forms of these frequencies as a point of 
comparison. We begin by we taking the exact expressions for E, L z , and Q given by Schmidt (Appendix B of Ref. 
f32|). expand to leading order in a, and then expand in 1/p. The result is 



E 



(1 



e 



2\2 



2p 



Li 



= m 2 M 2 p 



P 

s 2 (3 



8 M y/p 

2 ^ 2 - / 2(3 + e 2 ) 4(2 + e 2 )(3 + e 2 ) 



a 
M 



P 



3/2 



n5/2 



sm( 



(3.35) 
(3.36) 
(3.37) 



3 We have adjusted notation from Schmidt slightly to avoid confusion with our metric function 7. 
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The orbital frequencies become 



n r = 
n e = 
= 



,K 



K 



1 - 
1 + 
1 + 



3(1 -e 2 ) a 3(1 -e 2 ) sin* 



p M p 3 /2 

3e 2 a 3(1 + e 2 )sm9 min 
— ~ m pW 

3e 2 a 3(l + e 2 )sin6» min 

P 



M 



V 



,3/2 



+ 



a 2 



where 



is the Kepler frequency. 



,K 



1_ / 1-e 2 
M V p 



3/2 



(3.38) 
(3.39) 
(3.40) 

(3.41) 



C. Orbital frequencies of bumpy black holes 

Now examine how the orbital frequencies change if the black hole is bumpy. Begin with the Hamiltonian H. It 
remains conserved with value — m 2 /2, but its functional form is shifted: 



n 



1 



g af3 p a pp 



2 

H + Hi 



m 



(3.42) 



where Ti is the original (non-bumpy) Hamiltonian and TL\ gathers together the influence of the spacetime's bumpiness. 
To first order in b a p 1 



where 



Combining these expressions, we have 



which can be rewritten 



b afJ = g^g^b^ 



(3.43) 



(3.44) 



(3.45) 



Hi = --b^P v = -- 



dr dr 



(3.46) 



When we add bumps to a spacetime, shifting the Hamiltonian by Hi, the motion is no longer separable (except for 
the special case of equatorial motion) and the techniques used in Sec. IIII Bl for computing orbital frequencies do not 
work. However, since the spacetime is "close to" the exact black hole spacetime in a well-defined sense, the motion is 
likewise "close to" the intcgrablc motion. We can thus take advantage of canonical perturbation theory as described 
in, for example, Goldstein [41j, to calculate how the spacetime's bumps change the frequencies. 

The key result which we use is that the shift can be found by suitably averaging TL\ : 



m Slo 1 



dJ,, 



m ST 



(3.47) 



Notice that the derivatives are taken with respect to the action variables defined for the background motion; the 
averaging, which we denote with angular brackets, is likewise done with respect to orbits in the background. Once 
these derivatives are taken, it is simple to compute the changes to the observable frequencies. Expanding 



9J 



6uj 1 



ST 



(3.48) 
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we read out 



or 



5Q l 



sn i 



Suj 1 
t 



ST 



5u l 



r 2 

ST 
t 



(3.49) 



(3.50) 



The averaging used in Eq. (|3.47[) in described in detail in Appendix [AJ This procedure uses the fact that the radial 
and polar components of the background motion can be separated, and thus can be averaged independently. Following 
Appendix [XI this amounts to computing 



(Wi> = 



T*(27r) : 



dw r / dw 9 Hi [r(w r ), 6(w 9 )] T [r(w r ), 9(w 9 )] , 



(3.51) 



where T(r,6) is defined by Eq. (|3.4|) . T* is defined by Eq. (|Allj) . and where w r,e are angles associated with the 
separated r and 6 motions (defined and discussed in detail in Appendix [AJ . 

Before moving on, we discuss a few issues in practically computing (Hi) and the frequency shifts. Begin by 
expanding the Hamiltonian: 



Hi 



m 
~2 

~2 
1 

2E2 



dr 



ht\ -r-\ + Wr -r- 



dr 



d4> 



2b, 



dt dr 
dr dr 



2b, 



dtp dr 
dr dr 



2b. 



d<j) dt 
dr dr 



dr 



-r- I + bad, -r- 



d(f> dt 
£04>t — — 



[b tt T{r, 6) 2 + b rr R(r) + b ee Q{9) + b^(r, Of 



dr dr 
2b^(r,6)T(r,e)] 



(3.52) 



In going from the first line to the second, we use the fact that terms linear in dr/dr go to zero when we average since 
the radial motion switches sign after half a cycle. The final line of Eq. (|3 . 52[) is in a good form for averaging. 

In computing this average, we end up with (Hi) as a function of p, e, and # m i n . We likewise compute the actions 
J M using these parameters, and then compute the shifts to the frequencies and T using the chain rule. To set up this 
calculation, define an array bp which contains all the system's physical parameters: 



bp = (m,p,e,e min ) . 

Next, define the matrix J, the Jacobian of the actions with respect to these parameters: 

/ t n p _ dJ a 

W lot — a, 

obp 



Then, 



to* = , 
ob a 



ST 



d(H 



Ob, 



t 

a j 



(3.53) 
(3.54) 

(3.55) 
(3.56) 



where J 1 is the matrix inverse of the Jacobian J . 



IV. RESULTS I: ORBITS OF BUMPY SCHWARZSCHILD BLACK HOLES 



We now examine spacetimes and orbits for specific choices of ipi . Recalling that this function satisfies the Laplace 
equation, we take tpi to be a pure multipole in the 'Weyl sector" [i.e., in the coordinates of Eq. (|2.ip ]. As we will 
show in this section and the next, this smoothes out the bumps and cures the strong-field pathologies associated with 
orbits in the spacetimes developed in CH04. 

Note that a pure multipole in the Weyl sector will not correspond to a pure Geroch-Hansen moment of the black 
hole. For example, taking ip± to be proportional to an I = 2 spherical harmonic does not change only the moment 
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M 2 of the Geroch-Hansen sequence [Eq. (| 1 . 1 [) ] . However, it turns out that taking ip\ to be proportional to a spherical 
harmonic Y/o changes no Geroch-Hansen moments lower than M/ : taking ipi to be an I = 2 harmonic changes M 2 and 
higher moments; taking it to be an I = 3 harmonic changes M3 and higher; etc. A companion paper [28[ demonstrates 
this explicitly, and further shows that the dominant change for ipi oc Yiq is to the l-th Geroch-Hansen moment. 
Further, since the equations governing ipi and 71 are linear in these fields, one can choose ipi to be a combination of 
multipoles such that the resulting spacetime puts its "bump" into a single Geroch-Hansen moment. In this way, one 
can arbitrarily adjust the Geroch-Hansen moments of a spacetime (providing that the adjustments are small). 
We begin with spacetimes and orbits of bumpy Schwarzschild black holes. 



A. Quadrupole bumps (I = 2) 



First we examine an I = 2 perturbation in the Weyl sector. The perturbation -0i which satisfies Eq. (|2.2p and has 
an I = 2 spherical harmonic form is 



^ l r 2 (p,z) = B 2 M 



3 Y 2a (6 Wcyl ) = B 2 M 3 [5 3 cos 2 <9 W c y i 
(p 2 + z 2 ) 3 / 2 4 



- 1 



7T (p 2 + z 2)3/2 



(4.1) 



where cosf?wcyi = z/y'p 2 + z 2 . The dimensionless constant B 2 sets the magnitude of the spacetime's bumpiness 
for this multipole. Since we are treating the bumpiness as a perturbation, B 2 <gC 1. Transforming to Schwarzschild 
coordinates by Eqs. ([27]), (HU), (j2~T2)) . and (j2T3)l . we find 



ri= 2 (r,t 



B 2 M 3 [5 1 



7T d(r, 6) 3 



3(r - Mf cos 2 



d(r, ey 



- 1 



where 



d(r, 0) = (r 2 - 2Mr + M 2 cos 2 6) 1/2 



(4.2) 



(4.3) 



As a useful aside, the mapping from (p, z) to (r, 9) implies that any Weyl sector ipi can be transformed into 
Schwarzschild coordinates by putting 



2 1 2 
p + z 



cos 6 



Weyl 



d(r,0) , 
(r - M) 
d(r,6) 



cos 9 



Integrating the constraint (|2. 1 1|) and imposing the condition 71 (r — > 00) = gives 

(r-M) [c 20 (r) + c 22 (r) cos 2 6] 



1 [= 2 {r,6)=B, 



where 



C 2 oM 

C22 (r) 



d{r, Of 



2(r - Mf - 5M 2 (r - M) 2 + 3M 4 
5M 2 (r - M) 2 - 3M 4 . 



(4.4) 
(4.5) 



(4.6) 



(4.7) 
(4.8) 



Figure [T] shows the impact of an I = 2 bump on orbital frequencies as a function of p for a few choices of e and 
9 m . We show the three shifts 5Q r ' e '^, normalized by the bumpiness B 2 and rescaled by the asympotic weak-field 
dependence Sflf =2 oc p~ 7 / 2 (derived in Appendix iBj) . As we move to the weak field, the numerical results (solid lines) 
converge to the weak-field forms (dashed lines). The frequency shifts generically get substantially larger as we move 
into the strong field. The bumps have a very strong influence near the last stable orbit, plso = (6 + 2e), although 
the behavior is smooth and non-pathological. 



B. Octupole bumps (I — 3) 



Next, consider an I = 3 perturbation. In Weyl coordinates, we put 



^ 3 (p,,) = S 3 M^p^; (4.9) 
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= tt/3 



= tt/6 



= tt/3 







FIG. 1: Shifts to black hole orbital frequencies due to an I — 2 bump. The shifts <5JT' e ''* are normalized by the bumpiness 
parameter B2, and are scaled by p 7//2 ; this is because in the Newtonian limit, 5f2 r ' 8 '^ oc p~ 7 ^ 2 . The Newtonian limit (dashed 
lines) does a good describing the exact calculations (solid lines) for large p. This limit substantially underestimates the shifts in 
the strong field. Notice that the radial frequency shift changes sign in the strong field, typically at p ~ (10 — 13)M, depending 
slightly on parameters. This behavior is starkly different from the weak field limit. 



In Schwarzschild coordinates, this becomes 



B 3 M 4 1 



d{r, 9) 4 V tt 



5(r - M) 3 cos 3 6> 3(r-M)cos( 



d(r, ey 



d(r,9) 



From the constraint equation (12. lip and the condition 71 (r — > 00) — 0, we find 



B 3 M 5 17 



cw(r) + c 32 (r) cos 2 9 + c 34 (r) cos 4 9 + c 36 (r) cos 6 



where 



d(r, ey 



c 3 o(r) = -3r(r-2M) , 

c 32 (r) = 10r(r - 2M ) + 2M 2 

C3i(r) = -7r(r-2M) , 

c 36 (r) = -2M 2 . 



(4.10) 



(4.11) 



(4.12) 
(4.13) 
(4.14) 
(4.15) 



Notice that ipi and j[ =3 are proportional to cos 9. As such, their contribution to the averaged Hamiltonian (Hi) 
is zero: There is no secular shift to orbital frequencies from I = 3 bumps. This is identical to the result in Newtonian 
gravity, as discussed in Appendix |51 and holds for all odd values of I. 

As in the Newtonian limit, there will be non-secular shifts to the motion which cannot be described by our orbit- 
averaged approach. These shifts would be apparent in a direct (time-domain) evolution of the geodesies of spacetimes 
with odd / bumps. It would be a useful exercise to examine these effects and ascertain under which conditions odd I 
spacetime bumps could, in principle, have an observable impact. 



C. Hexadecapole bumps (I — 4) 



We conclude our discussion of Schwarzschild bumps with I — 4: 



, r 5 ^40(#Weyl) 



from which we obtain 



16 d(r,6») 5 V 7T 



( P 2 + z 2 y/ 2 ' 

35(r- Af) 4 cos 4 6» 30(r — M) 2 cos 2 ( 



d(r, 9) 4 



d(r, Of 



(4.16) 



(4.17) 
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Solving for 71 as before, we find 



where 



(r — M) C4 (r) + C4 2 (r) cos 2 9 + c^r) cos 4 



d(r, Of 



(4.18) 



C40M = 8(r - M) s - 36M 2 (r - M) b + 63M 4 (r - M) - 50M b {r - M) + 15M 8 , 
c 42 (r) = 36M 2 (r - M) 6 - 126M 4 (r - M) 4 + 120M 6 (r - M) 2 - 30M 8 , 
c 44 (r) = 63M 4 (r - M) 4 - 7(W 6 (r — M) 2 + 15M 8 . 



(4.19) 
(4.20) 
(4.21) 




FIG. 2: Shifts to black hole orbital frequencies due to an I — 4 bump. The shifts 5fl r ' are normalized by the bumpiness 
parameter B4, and are scaled by p 11 ' 2 , which sets the scaling in the Newtonian limit. As in the / = 2 case (Fig.[l}, exact results 
and the Newtonian limit coincide at large p, but there are significant differences in the strong field. The functional behavior of 
the radial frequency shift can be especially complicated in this case. 



Figure [5] presents the same orbits as are shown in Fig. [T] for / = 4. The shifts are normalized by the bumpiness 
-B4 and rescaled by the weak-field form 5Qf =i oc p~ n / 2 . The qualitative behavior is largely the same as for the 
quadrupole bump. In particular, we see once again that there are no strong-field pathologies in the orbit shifts, and 
that the degree of shift due to spacetime bumps is especially strong near the last stable orbit. The strong-field radial 
oscillations in 5tt r are even more pronounced than they were in the I — 2 case. This appears to be a robust signature 
of non-Kerr multipoles in the strong-field. 



V. RESULTS II: ORBITS OF BUMPY KERR BLACK HOLES 



We now repeat the above exercises on a Kerr black hole background. 



A. Quadrupole bumps (/ = 2) 



We begin with an I — 2 perturbation in the Weyl sector, starting with Eq. (|4.ip . Transforming to prolate spheroidal 
coordinates by Eqs. (|2.7[) and (|2.8|) , this becomes 



B 2 5 



3 cosh 2 u cos 2 v 



4 V Ti" \ sinh u sin v + cosh u cos 2 v 



— 1 I (sinh u sin v + cosh u cos v) 



-3/2 



(5.1) 



The corresponding 71 is 



"f[- 2 (u,v) = B 2 



cosh u [4 — cos 2v + (5 cos 2v — 1) cosh 2u + cosh 4m] 



(sinh 2 u sin 2 v + cosh 2 u cos 2 



5/2 



(5.2) 
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Following the logic of the Newman- Janis algorithm, we allow u to be complex, and replace cosh u with cosh u cosh u 
and sinh 2 u with (cosh u cosh u — 1). Making the coordinate transformation 



r — tacosa 
cosh u = — 1 , 



puts the result in Boyer-Lindquist coordinates: 



B 2 M 3 5 1 



= Box - 



4 V 7rd(r,(9,a) 3 
'5" r 



M 



3L(r, 9, a) 2 cos 2 



- 1 



d{r,9,a) 2 

L(r, 9, a) [c 2 o(V, a) + C22O*, a) cos 2 9 + c 2 i{r, a) cos 4 9] 
2 d{r : 9,a) 5 



where 



and 



d{r,9,a) = Vr 2 - 2Mr + (A-/ 2 + a 2 ) cos 2 6> , 
L(r,9,a) = y/(r - M) 2 + a 2 cos 2 9 , 



(5.3) 

(5.4) 
(5.5) 



(5.6) 
(5.7) 



(5.8) 
(5.9) 
(5.10) 



c 20 (r, a) = 2(r - Af) 4 - 5M 2 (r - A/) 2 + 3Af 4 , 
c 22 (r, a) = 5Af 2 (r - Af) 2 - 3A-f 4 + a 2 [4(r - A-/) 2 - 5M 2 ] 
c 24 (r,a) = a 2 (2a 2 + 5M 2 ) . 
Note that the result for can be found by taking the Weyl-sector perturbation given by Eq. (|4.ip and putting 

p 2 + z 2 -> d(r,9,a) , (5.11) 
i(r,9, (512) 



■ cos 9 



d(r, 9, a) 

This is the Kerr analog to the mapping described in Eqs. (|4.4|) and (|4.5 




FIG. 3: Shifts to Kerr black hole orbital frequencies for an I = 2 bump. As with the Schwarzschild results presented in Fig. 
[TJ the shifts SQ r ' e '^ are normalized by the bumpiness B2 and scaled by p 7 ^ 2 . Rather than examining a variety of orbital 
geometries, we here examine a few black hole spins, showing results for a — O.lAf, a = 0.5M, and a — 0.9AA Qualitatively, the 
results are very similar to what we find for the Schwarzschild case. The major difference is that the last stable orbit is located 
at smaller p, so that these orbits can get deeper into the strong field. The overall impact of the bumps is greater in these cases 
which reach deeper into the strong field. 



Figure shows how an I — 2 bump changes Kerr orbital frequencies. We focus here on how black hole spin affects 
our results, presenting results for a single orbit geometry (e = 0.3, 9 m = n/3). The major impact of black hole spin is 
to change the value of p at which orbits become unstable. For large spin, orbits can reach deeper into the strong field, 
accumulating larger anomalous shifts to their orbital frequencies. Aside from this behavior, spin has relatively little 
effect on the shifts: the three panels are similar to one another and to the Schwarzschild result (compare left-most 
panel of Fig. [TJ. Similar results hold for other orbital geometries, so we confine our plots to these results. 



17 



B. Octupole bumps (/ = 3) 



For the I = 3 Kerr bump, we begin with Eq. (|4.9[) . then follow the same procedure to take it to Boyer-Linquist 
coordinates as described for the / = 2 Kerr bump. The result is 



7i =3 M) 



i 



Trd(r,0,a) 4 



5L(r, 9, af cos 3 9 3L(r,9,a) cost 



d(r, 9, af 



BiM 5 pr 



C3o(r, a) + c 32 (r, a) cos 2 



d{r, 9, a) 
C34(r, a) cos 4 



c 36 (r, a) cos 6 



where 



d(r,6,a) 7 



c 30 (r,a) = -3r(r-2M) , 

c 3 2{r, a) = Wr(r - 2 M) + 2M 5 - 3a 2 

c 3 fr, a) = -7r(r - 2M) + 10a 2 , 

c 36 (r,a) = -2Af 2 - 7a 2 . 



(5.13) 
(5.14) 



(5.15) 
(5.16) 
(5.17) 
(5.18) 



As with the Schwarschild I — 3 bumps, V'i =3 an d 7i =3 are proportional to cos 9, so that (Hi) = 0. Thus, for Kerr as 
for Schwarzschild, there is no secular shift to orbital frequencies for / = 3, or any other odd value of I. We emphasize 
again that there will be non-secular shifts to the motion which our orbit-averaged approach misses by construction, 
and that it would be worthwhile to investigate their importance in future work. 



C. Hexadecapole bumps (/ = 4) 

For I = 4, we begin with Eq. (14.161) . Repeating our procedure to take this into a Kerr bump, we find 

3 , 



V4 =4 M) = 
7i =4 M) = 

where 

c 40 (r, a) 
c 42 (r,a) 

c 44 (r,a) 

c 46 (r,a) 
c 48 (r,a) 
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7T d(r, 9, a) E 



35L(r,0,a) 4 cos 4 (9 30L(r, 9, a) 2 cos 2 



d(r,9,a) 4 



d(r,0,a,y 



B A 



L(r, 9, a) [c 40 (r, a) + c 42 (r, a) cos 2 9 + c 44 (r, a) cos 4 9 + c 46 (r, a) cos 6 9 + c 48 (r, a) cos 8 9] 
8^ d(r,9,a) 9 



(5.19) 

- 1 
(5.20) 



8(r - M) 8 - 36Af 2 (r - Mf + 63M 4 (r - M) 4 - 50M 6 (r - M) 2 + 15M 8 , 
36M 2 (r - Mf - 126M 4 (r - M) 4 + 120Af 6 (r - M) 2 - 30A/ 8 
+2a 2 [I6(r - Mf ~ 54M 2 (r - M) 4 + 63Af 4 (r - M) 2 ~ 25Af 6 ] , 
63Af 4 (r - Mf - 70Af 6 (r - M) 2 + 15M 8 

+3a 2 [4Af 2 (9(r - A/) 4 - 21A/ 2 (r - M) 2 + 10M 4 ) + a 2 (l6(r - Mf - 36A/ 2 (r - Af) 2 + 21A/ 4 )] 



(5.21) 
(5.22) 



2a 2 [2a 4 (8(r - Mf - 9M 2 ) + 9a 2 M 2 (6(r - Af) 2 - 7Af 2 ) + 7A/ 4 (9(r - A//) 2 - 5Af 2 )] , 
a 4 (8a 4 



36a 2 A/ 2 + 63Af 4 ) 



(5.23) 
(5.24) 
(5.25) 



It is a straightforward exercise to numerically compute Sfl r ' e '^ using tp 1 ^ 4 and 7^ =4 and build the Kerr analogs to the 
results we show in Fig. [2l The results are not markedly different from those for Schwarzschild, modulo the fact that 
the orbit can (especially for large spin) reach deeper into the strong field and hence accumulate more bump-induced 
anomalous precession. Since there are no particularly surprising features in these results compared to what we have 
already shown, we will not show such plots. 



VI. SUMMARY AND FUTURE WORK 



This analysis significantly improves on the earlier presentation of bumpy black holes given in CH04, producing 
bumps with well-behaved strong-field structure, and extending the concept to Kerr black holes. These extensions 
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greatly expand the astrophysical relevance of these spacetimes. We have also demonstrated how Hamilton- Jacobi 
theory can be applied to orbits in bumpy spacetimes to categorize the anomalous precessions arising from their 
bumps in a reasonably straightforward manner. Bumpy black holes can now, at least in principle, be used as the 
foundation for strong-gravity tests with astrophysical data. 

It's worth re-emphasizing why we propose to use bumpy black holes, rather than using exact solutions which include 
black holes as a limit (for example, the Novikov-Manko spacetime [31| used in Ref. [30|). In large part, our choice is 
a matter of taste. Our goal is to tweak a black hole's moments in an arbitrary manner, so that the non-Kerr nature 
is entirely under our control. From the standpoint of formulating a null experiment, it arguably makes no difference 
whether the Kerr deviation takes one particular form or another. Any falsifiable non-Kerr form is good enough to 
formulate the test. To our minds, the nice feature of this approach is that if, for example, a theory of gravity specifies 
that a black hole should have the same moment structure as general relativity up to some I — L, it is simple to design 
a spacetime tailored to testing that theory. The case of Chern-Simons gravity discussed in Ref. [25[ shows that this 
motivation is not merely academic, but is motivated by plausible alternatives to general relativity. 

This analysis provides a complete description of "mass-type" bumps, i.e., perturbations to the mass moments Mi. 
We have not examined "current-type bumps," shifts to the spin moments Si. It is also worth bearing in mind that 
a pure multipole in the Weyl sector, %p\ tx Yiq(9) does not correspond to a shift in a pure Geroch-Hansen multipole. 
A companion paper by one of us (SJV) addresses these issues [281 ]. That analysis shows ho w sp in moments can be 
adjusted from their Kerr values by perturbing the Ernst potential describing the spacetime [27|. It also shows that 
a pure Weyl-sector multipole ip[ changes no Geroch-Hansen moments lower than /. As such, a pure Geroch-Hansen 
moment perturbation can be assembled by combining multiple Weyl-sector multipole perturbations. 

Aside from these more formal issues, the next major step in this program will be to use these foundations to 
formulate actual strong-field gravity tests that can be applied to astrophysical data. We imagine several directions 
that would be interesting to follow: 

• Extreme mass ratio inspiral (EMRI): The capture and inspiral of stellar mass compact into massive black holes 
at galaxy centers is one of the original motivations of this work. Much of the recent literature on testing and 
mapping black hole s pac etimes has centered on understanding the character of orbits in non-Kerr black hole 
candidate spacetimes [I5l[29ll43j. with an eye on application to gravitational- wave measurement of EMRIs. 

The full analysis of EMRIs in non-Kerr spacetimes is, in principle, quite complicated since their non-Kerr- 
ness breaks the Petrov Type-D character of Kerr black holes. As such, it may be quite difficult to accurately 
compute their radiation emission. It may not be quite so difficult in bumpy black hole spacetimes. Thanks to the 
smallness of their non-Kerr character, it may be fruitful to use a "hybrid" approach in which the short timescale 
motion is computed in the bumpy spacetime, but the radiation generation and backreaction is computed in the 
Kerr spacetime; a similar idea was suggested in the context of modeling EMRI events in Chern-Simons gravity 
[44| . Given that our goal is to formulate a null experiment, this hybrid may be good enough for a useful test, 
in lieu of solving the entire radiation reaction problem in non-Kerr spacetimes. 

• Black hole-pulsar systems: One of the goals of the planned Square Kilometer Array (SKA) |6| is the discovery 
of a black hole-pulsar binary system. If such a system is discovered, detailed observation over many years 
should be able to tease the multipole structure of the black hole from the data. Similar observations of neutron 
star-pulsar binary systems have already allowed us to make exquisite measurements of neutron star properties 
and gravitational- wave emission [45j . The tools developed here may already be adequate for doing this analysis 
since such binaries will have a relatively slow inspiral time. 

• Accretion flows on black hole candidate: Programs to observe the (presumed) black hole at the center of our 
galaxy are maturing very quickly; programs to study accretion flows onto stellar mass black holes in x-rays are 
already quite mature. In our galactic center, the most precise measurements come from millimeter wavelength 
radio emission from gas accreting onto this central object. The precision of these measurements is increasing 
to the point where we will soon be able to use them to map the detailed strong-field spacetime structure of the 
spacetime near Sagittarius A*. It would be a worthwhile exercise to repeat analyses of the appearance of these 
flows in Kerr spacetimes [46j for bumpy spacetimes to see how accurately such measurements may be able to 
probe the central object's multipoles. Such an analysis will require developing imaging maps of bumpy black 
holes, going beyond the orbit frequency analysis we have given here. In the domain of x-rays, it would be very 
interesting to extend work on quiescent accretion flows (as in, for example, Ref. [5j) to include accretion in 
bumpy black hole spacetimes. It may be possible to extend those analyses to fit for multipoles beyond the black 
hole's mass and spin. 
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APPENDIX A: AVERAGING FUNCTIONS ALONG BLACK HOLE ORBITS 



Key to computing the shifts that a black hole's bumps impart to its orbital frequencies is averaging the Hamiltonian 
perturbation T~L\ = Hi(r, 9) along the orbit. The averaging we wish to use is a time average: 

(Wi)= lim -L f Hi[r{t),B(t)\ dt . (Al) 

Implementing this integral is somewhat involved, so we review the procedure here. These details were first developed 
by Drasco and Hughes (Ref. [loj]), hereafter DH04. 

Most GR textbooks (e.g., Ref. @, Chap. 33) describe the motion of a small test body orbiting a Kerr black hole 
using Eqs. (|3.1[) - (|3 .4|) . Although the coordinate motions are formally separated in these equations, in a practical 
application they are not quite "separated enough" thanks to the factors of £ = r 2 + a 2 cos 2 9 which appear on the 
left-hand side of these equations. These factors couple the radial and polar motion, and complicate averaging functions 
of the form /(r, 9) that are computed on an orbit. Unless the r and 9 periods coincide (or are in an integer ratio) one 
cannot easily average over the r and 9 motion. 

This residual coupling is eliminated by parameterizing the orbits using what is now often called "Mino time" A, 
defined by dX = dr/S. Equations (|3.1[) - (|3.4[) become 

d r y 

d9^ 2 



R(r), 


(A2) 


0(0), 


(A3) 


$(r,0) , 


(A4) 


T(r,9) . 


(A5) 



d\J 
(W_ 

rfA 

d<t> 

dX 
dt 

dX 

Since this time variable explicitly separates the r and 9 motion, it is simple to construct r(X) using Eq. (|A2[) . and 
likewise to construct 9(X) with (|A3j) . It is also simple to compute the r and 9 periods in Mino time: By inspection of 
Eqs. (|A2[) and (|A3|) . and taking into account symmetries of the motion and the appropriate functions, we have 

F" dr 

A r = 2 / -f = , (A6) 
A e = 4 -¥= (A7) 



Following DH04, we define frequencies conjugate to these periods, 

T r ' e = 2ir/A r ' e , (A8) 

and then introduce angles 

w r,9 = T r,9 x ( A9 ) 
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We then write the radial motion as a function of w r , 

r(w r ) = r(A = w r /T r ) ; (A10) 

we likewise parameterize the polar motion using w e . The key concept behind the averaging is that we now allow the 
angles w r and w 6 to separately vary. This allows us to separately average the r and 9 motions. 

We are nearly ready to use these tools to average our Hamiltonian. Before doing so, we define T', the time- function 
T(r, 9) averaged over the angles w r and w e : 

T'e— -/ dw r / dw e T[r(w r ),9(w 9 )] . (All) 
(2tt) Jo Jo 

[In DH04, this quantity was denoted by a capital gamma; we are adjusting the notation slightly to avoid conflict with 
r as defined in Eq. (|3.34p . It is something of an abuse of our notation, since the time t is not periodic and T l typically 
denotes a frequency] 

As shown in DH04, it is now simple to compute the long-time average of any black hole orbit functional: 
(/) = km -L f f[ r (t),e(t)]dt 

T— >oo Zl J _rp 

dw r dw e f[r(w r ),9(w e )]T[r(w r ),9(w e )]. (A12) 



TW Jo 

This is the procedure we use to compute (Hi) for all the computations presented in this paper. 

APPENDIX B: NEWTONIAN PRECESSION FREQUENCIES 

For weak-field orbits, we expect that the bumpy black hole frequency shifts [Eqs. (|3.47|) - (|3.50|) ] are well described 
using Newtonian gravity. In this appendix, we compute the relevant Newtonian frequency shifts; in Sees. HVl and lv | 
we show that our general formulas limit to the results we develop here for large radius orbits. 

As in our relativistic calculation, we compute frequency shifts due to multipolar "bumps" by examining the variation 
of a perturbed Hamiltonian with respect to an orbit's action variables: 



mmi= mxi (B1) 



(Since there is no distinction between coordinate and proper time in Newtonian gravity, these Sfl 1 are the measurable 
frequencies.) The actions arc defined just as in the relativistic case, Eqs. (|3.14p - ({3.160 . For a body of mass m 
orbiting a mass M in Newtonian gravity, they become 



J r = m ^p^U—L^ -lj , (B2) 

Je = m^fp~M (1 - sin 9 m ) , (B3) 

J0 = m^ppM sin 9 m . (B4) 

The perturbation to the Hamiltonian for orbits is proportional to the perturbation to the potential: 

Hf" wt = mSVt (r, 9) = "^f^ o (cos 9) . (B5) 

To perform the averaging, we first reparameterize both the radial and angular motion in a manner similar to what 
we use for black hole orbits: 

r=r^— , (B6) 
1 + e cos ip r 

cos 9 = cos9 m cos (tp r — xo) • (B7) 
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Notice that the radial and angular motions vary in phase with one another in the Newtonian limit: both complete a 
full cycle as the angle ip r varies from to 2tt. The angle Xo ls an offset phase between these motions. Inserting this 
parameterization into the Newtonian equations of orbital motion, we find 




ecos ip r y 



(B8) 



The averaged Hamiltonian is then given by 



( ^ CWt> = TWo 6Vl[r ^' m]dt 



rn 



2 7T 



dlp r 

~~dt 



5Vi [r(^v),cos0(Vv)] dip r 



(B9) 



where Tk = 2ttM^p 3 /(1 — e 2 ) 3 is the Keplerian orbital period. [To put this in a more familiar form, recall that the 
orbit's semi-major axis A = pM/(l — e 2 ).] 

Using these results, we now examine the Newtonian limit of I — 2, 3, and 4 black hole bumps. 



1. Quadrupole (I — 2) 

The quadrupole bump is given by the potential 



for which we find 



(Hi 



SV l = 2 = 



B 2 M 3 
4r 3 



cos 



1] 



1 



ZB 2 


1 


8M 


p 7/2 


3B 2 


1 


8M 


p 7/2 


W 2 


1 


8M 


p 7/2 



-^(l-e 2 ) 2 (3sin 2 9 n 



Varying this averaged Hamiltonian with respect to J r 
Sil r = 

sn e = 

(50^ = 



■ e ) [1 — 3sm 
we find 



b (l_ e 2)3/2 
7T 



sin 2 Or, 



(1 



,2^,3/2 



sin 



3V1 



3V 7 ! - e 2 



2sin( 



Vl-e 2 - 1 



(BIO) 
(Bll) 

(B12) 
(B13) 
(B14) 



These frequencies are written using the Kepler frequency lok — 2tt/Tk- 

Equations (|B12|) - (|B14|) reproduce well-known results for motion in a spherical potential augmented by a quadrupole 
perturbation. To facilitate comparison with the literature, it is useful to change our description of the orientation of 
the orbital plane from 9 m (the minimum angle 9 reached over an orbit) to i = 7r/2 — 9 m (the orbit's inclination with 
respect to the equatorial plane). We then construct the precession frequencies 



n aps 



sn e - SQ r 



to K x 



n planc = (5ft* - sn e = -v" x 



,K 



3B 2 1 

8M p 2 

3B 2 1 

4Af p 1 



(5 cos 2 i - l) 



5 

— COS L 
7T 



(B15) 
(B16) 



f^apsis describes the frequency of the orbit's apsidal precession within its orbital plane, and O plane the frequency at 
which the orbital plane precesses around the symmetry axis. These frequencies reproduce expressions that can be 
found in the literature; cf. Sec. 12. 3C of Ref. [4l|, and Sec. 3.6.2 of Ref. |42j |. 

The apsidal precession (|B15|) is due to beating between the orbit's radial and polar motions. In many relativity 
applications, one is interested in the precession of an orbit's periastron, which is due to a beat between the orbit's 
radial and azimuthal motions. For the quadrupole bump, we find 



sn<t> - sn r = uj k x 



3B 2 1 
8M p 



(5 cos 2 l — 2 cos u — l) 



(B17) 
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Taking the equatorial limit (t = 0), we see that this result agrees with Eq. (A4) of CH04 provided we identify their 
parameter Q with B2M 3 ^/5/4ir. Comparison of the potential used in CH04 [their Eq. (A2)] with our form [Eq. (|B10[) ] 
shows that this identification is exactly correct. 



2. Octupole shift (I = 3) 

The octupole bump is given by the potential 



SV 



1=3 



B 3 M 4 17 
4r 4 



[5 cos 3 0- 3 cos 6} , 



leading to 



{Hl) = \fl (1 " ^ COS 9m (5 C ° S2 9m ~ 4) C ° S X0 ■ 



(B18) 



(B19) 



This is proportional to cosxoj and so depends on the phase offset of the radial and angular motions. Over very long 
timescales, precessions will average out this dependence. However, on timescales that are not long enough for xo to 
vary over its full range, there is a residual impact, leading to octupolar shifts to observable frequencies. This can 
be very important in Newtonian celestial mechanics. For black hole applications, where the r and motions do not 
vary in phase with one another, this averaging will be much stronger, and treating (H\) — should be much more 
accurate. 

Varying the averaged octupole Hamiltonian shift, wc find 



sn r = 

SQ e = 

sn e = 



3£ 3 1 7 



16eM p 9 / 2 V 7T 
■Mh 1 / 7 



-(1 - e 2 ) 2 (l - 4e 2 ) (5 cos 2 9 m - 4) cos X o , 



16eM p 9 / 2 V 7T 
— 4e 2 sin 9 m tan 



(1 - e 2 ) 3 / 2 



cos 6 



4 - v 7 ! - e 2 + 4e 2 (l + ^fl-e 2 ^ + 15e 2 si 
\]\ -e 2 + 4e 2 + v 7 ! - e 2 ) 



sin 2 Or, 



COS Or, 



cosxo , 



3^3 1 /7^ 1 _ e 2j 3 /a 



16eM p 9 / 2 V 7T 
-4e 2 (sin6' m - 1) tan# 



cos 



4 - 4 VI - e 2 + 16eVl - e 2 + e 2 [16 + 15 (sin6» m - 1)] sinf 



1 - y/l - e 2 + 4e 2 f 1 + v 7 !^ 



cos 0„ 



cos xo 



(B20) 



(B21) 



(B22) 



As expected from the form of the averaged Hamiltonian, all precession frequencies are proportional to the cosine 
of the offset phase xo- On average we therefore find no influence from the I = 3 perturbation, nor from any odd I 
multipolar bump. The "instantaneous" impact of odd I bumps is, however, non-zero. 



3. Hexadecapole shift (I = 4) 

Finally, for the hexadecapole bump, we have 

H=4 _B 4 M 5 [9 



SV 



16r 5 



[9 

■J- [35 cos 4 6»- 30 cos 2 6> + 3] , 



(B23) 



leading to 



(Hi) = ^-4\ - (1 - e 2 ) 3/2 [8(2 + 3e 2 ) - 20 cos 2 6» m (4 + 6e 2 + 3e 2 cos 2xo) + 35 cos 4 6» m (2 + 3e 2 + 2e 2 cos 2xo)l • 
256p° V 7r 

(B24) 

We focus on the secular (long-time average) precessions, and average this over xo, leaving 



3mB 4 19 
256p 5 



(1 - e 2 ) 3/2 (2 + 3e 2 ) (8 - 40 cos 2 m + 35 cos 4 0„ 



(B25) 
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The precession frequencies which arise from this are 



-S^v! (1 " e2)2(8 " 4Ocos20m+35cos40m) ' (B26) 

155 4 1 



^(1 - e 2 ) 3/2 [8 (8 + 3e 2 (3 + VT^)) - 4 (62 + e 2 (63 + 30^1^)) cos 2 ft, 



256Mp n /2 V 7T 
+7 (28 + 3e 2 (9 + 5^1 - e 2 ) ) cos 4 ft 



= ^ + 2S^/f (1 - e2)3/2(2 + 3e2)sin ^( 28 - 40C ° s2 ^ ' 



(B27) 
(B28) 
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